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Abstract 

The effect of multibody massless neutrino exchanges between neutrons 
inside a finite-size neutron star is studied. We use an effective Lagrangian, 
which incorporates the effect of the neutrons on the neutrinos. Following 
Schwinger, it is shown that the total interaction energy density is computed 
by comparing the zero point energy of the neutrino sea with and without 
the star. It has already been shown that in an infinite-size star the total 
energy due to neutrino exchange vanishes exactly. The opposite claim that 
massless neutrino exchange would produce a huge energy is due to an improper 
summation of an infrared-divergent quantity. The same vanishing of the total 
energy has been proved exactly in the case of a finite star in a one-dimensional 
toy model. Here we study the three-dimensional case. We first consider the 
effect of a sharp star border, assumed to be a plane. We find that there is a 
non-vanishing of the zero point energy density difference between the inside 
and the outside due to the refraction index at the border and the consequent 
non-penetrating waves. An analytical and numerical calculation for the case 
of a spherical star with a sharp border confirms that the preceding border 
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effect is the dominant one. The total result is shown to be infrared-safe, thus 
confirming that there is no need to assume a neutrino mass. The ultraviolet 
cut-offs, which correspond in some sense to the matching of the effective 
theory with the exact one, are discussed. Finally the energy due to long 
distance neutrino exchange is of the order of 10 -8 -10 -13 GeV per neutron, 
i.e. negligible with respect to the neutron mass density. 
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I. INTRODUCTION 



The massless neutrino exchange interaction between neutrons, protons, etc., is a long- 
range force [[[]- 0. In a previous work Q, the long-range interaction effects on the stability 
of a neutron star due to multibody exchange of massless neutrinos have been studied. We 
have shown that the total effect of the many-body forces of this type results in an infrared 
well-behaved contribution to the energy density of the star and that it is negligible with 
respect to the star mass density. This is in agreement with two recent non-perturbative 
calculations done by Kachelriess || and by Kiers and Tytgat 0. 

This work is in contradiction with the repeated claim by Fischbach J7| that, unless the 
neutrino is massive, neutrino exchange renders a neutron star unstable, as the induced self- 
energy exceeds the mass of the star because of the infrared effects associated to neutrino 
exchange between four or more neutrons. In our opinion, the latter "catastrophic" result 
is a consequence of summing up large infrared terms in perturbation outside the radius of 
convergence of the perturbative series. The non-perturbative use of an effective Lagrangian 
immediately gives the result without recourse to the perturbative series, and the result is 
small. 

Smirnov and Vissani ||, following the same method as in ref. 0, summing up multibody 
exchange contributions order by order, showed that the 2-body contribution is damped by 
the blocking effect of the neutrino sea ||. They guessed that this damping would apply 
to many-body contributions, and hence would reduce the catastrophic effect claimed by 
Fischbach. In our previous work |4| , we also considered the effect of the neutrino sea inside 
the neutron star. This effect has been introduced in our non-perturbative calculation by 
using Feynman propagators of neutrinos inside a dense medium, which incorporate the 
condensate term. We noticed that this condensate is present, but in our opinion it is not the 
most important of the effects that neutralize the catastrophic effect expected by Fischbach, 
since it only brings a tiny change to the non-perturbatively summed interaction energy, 
leaving unchanged our conclusion that the weak self-energy is infrared-safe. 

We have also stressed in J| that the neutrino condensate was related to the existence of 
a border Q This was demonstrated in the (1 + 1) -dimensional star in [II]]: the blocking 
effect, which implies the trapping of the neutrinos inside the star while the antineutrinos are 
repelled from it, is a natural consequence of the existence of a border. Indeed a proper treat- 



At that point, we should call the reader's attention to a minor mistake that was made in the 
previous calculation Q: a pole was forgotten in the calculation of the weak self-energy, and its 
contribution is exactly annihilated by the condensate's, as shown in ref. |i0| : the computation of 
the self-energy gives zero when the forgotten pole is taken into account. 
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ment of the effect of the border automatically incorporates the condensate contribution as a 
consequence of the appropriate boundary conditions for the neutrino Feynman propagator 
inside the star. 

Since our treatment directly incorporates the neutrino sea effect and, on the other hand, 
since we stick to our strategy of directly computing the total neutrino interaction energy 
by using an effective Lagrangian, we have in a sense thus generalized the result of Smirnov 
and Vissani || because our result holds in a non-perturbative way and accounts for all the 
n-body contributions, while theirs holds for 2-body contributions only. 

It has been objected JTTJ] to our study 0] that we worked in an approximation where we 
neglected the border of the neutron star. Our belief is that this simplifying hypothesis does 
not change the fundamental result that the total effect of the multiple neutrino exchange to 
the energy density of the star is not infrared- divergent. 

Indeed, in []10] we have proved that, in (1 + 1) dimensions with borders, the result of 
|D for the infinite star without borders is kept unchanged: the net interaction energy due 
to long-range neutrino exchange is exactly zero. This result was obtained by computing 
Feynman vacuum loops with neutrino propagators derived from the effective Lagrangian, 
which incorporates the neutron interaction. The latter propagators are not translational- 
invariant, because of the star borders. We have found a physically simple explanation for 
the vanishing of the net interaction energy. It relies on the fact that the negative energy 
states, in the presence of the star and without the star, are in a one-to-one correspondence 
and have exactly the same energy density. It results that the zero point energy is the same 
with and without the star, for any density profile of the star. 

The main goal of the present work is to follow on taking into account the finite-size and 
border effects. Two main conclusions of ref. |T0| are useful for the (3 + l)-dimensional star: 
i) the natural connection between the neutrino sea and the border, ii) the correct definition 
of the zero-energy level of the Dirac sea. 



From |TIJ, we know that the zero-energy level of the Dirac sea has to be adjusted by 
comparing the asymptotic behaviour of the wave functions far outside the star with the 
free solutions in the absence of the star. From there, we know the correct ie prescription, 
which has to be imposed in the propagators of the neutrinos in the presence of the star; we 
could in principle compute the closed loops to get the vacuum energy density. However, this 
approach is technically very difficult. A simpler method, the derivation of which is recalled 
in section [TI], is to simply add the energy density of the negative energy solutions in the 
presence of the star, minus the same in the absence of the star. 

The vanishing of the neutrino exchange energy density in the star, found in the case of 
an infinite star M and in (1 + 1) dimensions [I0|, which will be summarized in section [Til], is 
not valid in (3 + 1) dimensions. The main reason for that will be illustrated in section |TV|by 
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zooming to the border effect, i.e. considering a plane border in (3 + 1) dimensions. There 
is a non-trivial refraction index that modifies the wave energy densities as they penetrate 
the star. Some waves are forbidden to penetrate and this yelds the dominant contribution. 

In section [V], we perform analytical and numerical calculations, which take into ac- 
count the curvature of the border and use a spherical star with a sharp border. We find 
a very simple approximate formula for the zero point neutrino energy density in the star, 
and demonstrate numerically the validity of this approximation. It results that indeed the 
neutrino-induced energy density in the star does not vanish and is dominantly explained by 
the above-mentioned border effect. 

In any case, these non-vanishing neutrino exchange energy densities are all perfectly 
regular in the infrared and do not present any resemblance to Fischbach's effect. On the 
other hand they are ultraviolet-singular. This is not unexpected since anyhow our effective 
Lagrangian is only valid below some energy scale where the neutrons may be considered at 
rest. In section |VI], we also discuss in some detail the effect of decoherence for distances 
larger than the neutrino mean-free path, which smoothes down the ultraviolet singularity. 

We conclude that the stability of compact and dense objects such as a white dwarfs, 
neutron stars, etc., are not affected by the neutrino exchange, even if neutrinos are massless. 



II. GENERAL FORMALISM 

In our calculations, we assume that the material of the neutron star is made exclusively of 
neutrons, among which neutrinos are exchanged. The density-dependent corrections to the 
neutrino self-energy result, at leading order, from the evaluation of Z°-exchange diagrams 
between the neutrino and the neutrons in the medium, with the Z° propagator evaluated at 
zero momentum. The vacuum energy-momentum relation for massless fermions, E = \q\, 
where E is the energy and \q\ the magnitude of the momentum vector, does not hold in 
a medium In our case, following refs. |13| and [Kj, they can be summarized by the 



following dispersion relations: 



E u = q = \q\+b for \q\ > \b\ 

V\q\ or, (1) 
for \q\ < \b\ 



Eh, 




where 

b ~ -V2G F n n /2 ~ -0.2 10~ 7 GeV ~ -lO^fm -1 , and n n ~ 0.4 far 3 . (2) 
In this paper we will use a star radius of 

R ~ 10 km, whence bR ~ -10 12 . (3) 
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In eq. flU), b summarizes the zero-momentum transfer interaction of a massless neutrino with 
any number of neutrons present in the media. Sensibly enough, it depends on the neutron 
density which will be assumed to be constant for simplicity. When a neutrino sea is present 
H Q, the neutrino condensate does not sensibly modify the value of b ||. 

In the approximation where the neutrons are static, in the sense that they do not feel 
the recoil from the scattering of the neutrinos, we can study the neutrino exchange through 
the following effective Lagrangian as done in ref. [|J: 

£ eff = iV$ u(r) - b 77 7o z/ Q[R - r) , (4) 

where R is the radius of the neutron star. 

The dispersion relations in eq. ([!]) show a displacement of the energy levels for the 
different modes, a negative shift for neutrinos, and a positive one for antineutrinos; the 
Dirac sea level is displaced. Would the neutron star occupy the whole Universe, it would 
just mean a change of variables, with no physical consequence. The finite size of the star 
changes the picture. Notice that b acts as the depth of a potential well, as it will be 
considered in section |V|. It is repulsive for antineutrinos and attractive for neutrinos, which 
then condense. 

A traditional way to estimate the energy induced by neutrino exchange is to compute 
first the exchange potentials involving 2, 3, 4, ... neutrons and then add their contributions 
integrated over the neutron positions in the star. This is the way chosen by Fishbach et al. 
to compute the weak self-energy. The drawbacks of this method are the following: 

i) The calculation turns out to be so difficult that many approximations are necessary. 

ii) It assumes implicitly that every neutron may interact only once with neutrinos, iii) The 
resulting interaction energy increases with the number n of neutrons involved grossly as 
(bR) n ~ 10 12n and eventually becomes very large. In fact it would go to infinity if these 
authors did not stop when n equals the number of neutrons of the star, as a consequence of 
the above-mentioned unjustified assumption that a given neutron cannot interact more than 
once. The large parameter bR (eq. @) reflects the infrared-sensible behaviour of each term 
in the series. The main problem of this approach is that the summation is done outside the 



2 b is given by: 

i— , . , 8\/2G\F/i . . 
b ~ -V2G F (n n - n v )/2 — — n v {E v ) , 

where {E v ) is the average neutrino energy of the medium, defined in its rest frame, and n v the 
neutrino density {n v ~ 4 x 10 -23 fm -3 ). These corrections to the fermion propagation would give 
higher-order effects, and we disregard them in the present work. 
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radius of convergence of the perturbative series, leading to an unacceptably huge result J7j. 
In ref . 0] , we used instead a simpler and more direct method, which is non-perturbative and 
based on an effective action. This method does not involve uncontrollable approximations, 
it does incorporate automatically multiple interactions of one same neutron and finally it 
leads to a totally reasonable result. This method is also followed in recent works, refs. || 
and 0. We should insist that, in this approach, there is no extra assumption added to the 
problem. 

In our work, we use the Schwinger tools JT4J in order to compute the density of weak 



interaction energy w(x) due to the multibody neutrino exchange. It is given by the difference 
between the energy density for a neutrino propagating in the "vacuum" defined by the 
neutron star, |6), and the corresponding one for the real, matter-free, vacuum, |0): 

w(x) = (6\H(x)\6) - (0\H (x)\0) , (5) 

where 7i{x) and 7io{x) are the Hamiltonian densities of free and interacting neutrinos. 
Concretely, to compute this weak interaction energy analytically, using , we can write it 
as: 

d 

w(x) = -^Tr-^tr {70 [S F {x, y) - S { p\x, y) } . (6) 



The r.h.s. does not depend on x° from time-translational invariance. 

In diagrammatic formulation, the formal eq. corresponds to the computation of the 
difference of the diagrams in fig. [l|. 




q Yo Wo 

FIG. 1. Schematic representation of w, eq. ^). The simple and double line represent the 
neutrino propagator, in the matter-free vacuum and in the vacuum of the interacting medium, 
respectively. 

Using the notation 
and from the definition 



t 



E 

jz,l,» 



dE, 
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S F (x, y) l0 = 9(x° - y°) t, ^ n+ (xM + (y) e~ iE ^^ 

- %° - x °) % *«-(*)*L(fl e-^-^-^ , (?) 

where E n _ (E n+ ) is the negative (positive) energy of the eigenstate: H^/ n± = E n± ^/ n± , w(x) 
can be rewritten fl4j as: 



w\x) 



t E n _¥ n _(x)* n _(x) - t E^\x)¥°!(x) , (8) 

J n— Jri- 



where (0) refers to the matter-free vacuum and where we have taken the limit y° —>■ x° with 
y° > x°. Taking the average of the limits y° — > x° with y° > x° and y° — > a; with y° < x°, 
as done in ||, would lead to: 

w sym (?) = \$: E n ^_(x)* n _(x)-~^ E^°\x)¥°!(x) 

~ \ t E n+ (x ) * n+ (x ) + It, <> < 0) (x) *«g (5) • (9) 

The states are normalized in such a way that the two densities of states p{E), with and 
without the star (free case), coincide in the asymptotic region \x\ — > oo so that w(\x\ — > 
oo) = 0. 

For later use let us just remind the reader that the same result (|8]) may be obtained from 
the time Fourier-transformed propagator: 

S F (x,y,q ho = ^ £ g o _ K+ + , £ + g o _ En _ _ Z£ • ( 10 ) 

The result (H) is obtained from fllQl) by integrating — q° Sf{x, y, g°)/(27r) on the complex 
q° plane, by closing the contour on the upper half. Obviously an appropriate choice for the 
ie convention is crucial here, as discussed at length in [[K]]. It is also clear that formula (§) 
leads to much simpler calculations than the direct calculation of the loop in (|6]). 

Obtaining the weak self-energy of the finite neutron star is equivalent to calculating the 
neutrino propagation in a background of neutron fields density with a border. We did that 
analytically in (1 + 1) dimensions with a sharp border, and in (3 + 1) dimensions with a 
flat border. The spherical symmetry has been studied both analytically and numerically. 



III. THE TOY EXAMPLE OF (1 + 1) DIMENSIONS 

Working directly with eq. (|B|) requires an interacting Feynman propagator that takes 
the existence of a border into account. Doing this in (3 + 1) dimensions is a big task. 

It is feasible and theoretically fruitful to perform the analytical calculation of eq. (|6]) in 
(1 + 1) dimensions as a toy model. The knowledge extracted from that study will help us 
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face the (3 + 1) realistic case by using the general formalism anticipated in section |T[ We 
will now summarize the (1 + l)-dimensional results presented in more detail in [ [lOfl . 

The (1 + l)-dimensional toy model is presented in ref. |1U|| . We will summarize the 
computation and the result in order to extract useful information for the more realistic case. 

We consider (1 + 1) massless fermion Feynman propagators in a space with two regions 
separated by a border. The two regions, inside and outside the star, have different fermion 
dispersion relations, as seen in eq. ([[]). 

For simplicity, we consider only one sharp border located at z = and use an effective 



neutrino Lagrangian, which summarizes the interaction with the neutrons 15 



C cS = iV{z)pv{z)-b9{z)v{z) 1 %{z) , (11) 

where b is given by eq. @. The details of the computation of the interacting propagator 
can be found in ref. [ITH]. In momentum space, the propagator can be written as follows: 



S F (qf, q l ) = 2tt5 (q{ - ql) -4- + f -j-A 



where: 



^ = ? = ^ T -With ( q< T = q» = (q ,q) , 

^ = ^-with ( q> r=(q -b,q) ; (13) 

with s = sign(g )sign(g — b), and a z = 7°7 1 . The sign s results from the appropriate 
Feynman boundary conditions and from the right choice of the zero energy level. It gives an 



adequate time convention. We should emphasize that the propagator ( |T2| ) is infrared-safe, all 
the jf s in the denominators being regularized by Feynman's prescription for the distribution 
of the propagator poles in the complex-g plane |L0] . 



The expression for l/(^)*, given by eq. ([13]), can be appropriately rewritten as 

' 1 



— + 2mi>8 Iqi) 6(-qo)6(q - b) . (14) 



In ref. [|l(| , we have demonstrated that the main contribution to the weak self-energy comes 
from the first term of the propagator flT2p. We thus identified the expression given in 
eq. (|i~4D with the one for the effective propagator, which takes into account the "bulk" of the 
neutron star. As a matter of fact, the effective propagator for the infinite star [f|] coincides 
with the propagator (|H|), except for the time convention introduced by the proper boundary 
conditions, responsible for the second term in the l.h.s. of eq. (|1~4|). This term is nothing 
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else than the condensate contribution, i.e. the Pauli blocking effect of the neutrinos trapped 
into the star by the attractive potential. In refs. and f8|, the same condensate term in 
the l.h.s. was introduced by hand. Equations ( |TB| ) and (|T^) give a confirmation of the idea 
proposed in refs. [|J and [|15j: the condensate is a consequence of the existence of a border. 

The condensate is physically understandable. As we have tuned the level of the Dirac sea 
(see eq. ([!])) outside the star (to the left), and as our states extend over all space, far inside 
the star (to the right), the level corresponds to filling a Fermi sea above the bottom of the 
potential b < q < 0. This obviously induces a Pauli blocking effect and eq. ([!]) anticipates 
this result: \b\ is a lower bound for the momentum of the positive-go states. 

Now, with the computation of the propagator, it is easy to calculate the weak self-energy 
density w(z) for this (1 + l)-dimensional star by following eq. @. Some details about 
this rather cumbersome calculation may be found in pTJfl . It is worth while emphasizing 



that the vanishing result (^) for the border contribution requires the change in the order of 
integration in the momentum space over the variables qo, q\ and q{. The latter is only possible 
in the framework of an implicit ultraviolet regularization scheme. As a matter of fact, the 
choice of a certain order of integration over the variables is an ultraviolet regularization 
method. By integrating over the momenta we found that this difference vanishes exactly: 

w(z) = . (15) 

This vanishing interaction energy in (1 + 1) dimensions has a simple physical explanation: 
the presence of the border does not disturb the wave functions (up to a phase). There is a 
one-to-one correspondence of the negative-energy states in the two media z < and z > 0, 
and by interchanging the sum and integral as an ultraviolet regularization procedure, each 
term in eq. (|8|) vanishes exactly. In other words, massless propagating neutrinos are not 
reflected by the border of a (1 + 1) star, and the probability density is not affected by the 
border either. 

The use of eq. (§) allows a trivial generalization of the vanishing result flip] ) in (1 + 
1) dimensions to more complicated structures, for example a star with two sharp borders 
or, still better, to a continuously varying neutron density, as long as the density vanishes 
asymptotically outside the star. 



IV. FLAT BORDER IN (3 + 1) DIMENSIONS 

The computation of the neutrino loops (fig. [I]) in the (3 + l)-dimensional case, with a 
finite star, as was done in the preceding section for (1 + 1) dimensions, is technically rather 
difficult, if only because the neutrino propagators are not simple. In the preceding section 
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we have seen that the use of eq. @ simplified the calculation a lot, reducing it practically 
to triviality in the latter case. We will therefore stick to it in the (3 + 1) case. 

We might have faced the (3 + l)-dimensional problem analytically with interacting 
propagators and the approach that led to eq. (]6|), by using a simplifying assumption: the 
influence of the border can be reasonably neglected if the weak energy density is computed 
far inside the star. This amounts to considering a star large enough for the contribution 
of the bulk effective propagator ( p!4|) to be the main one. In other words, it would amount 
to assuming an infinite star but taking into account the existence of a matter-free vacuum 
infinitely far from the star centre. This matter-free vacuum allows us to fix the zero energy 
level, and consequently the adequate time convention for the effective propagator Jl^ . 

The latter calculation for the infinite star has indeed already been performed in our 
previous work ||. There, w(x) was computed using the Schwinger-Dyson expansion of 
eq. (16) in |4[] and the Pauli-Villars procedure as a means of regularization. As already 
noted in the introduction and in ref. fTPf , a pole (go = \o\ + b, in the case \q\ < b) has been 
forgotten (see fig. 1 in |TIJ) in the analytic continuation, which led to a wrong non-null result 
for w(x) when the neutrino condensate effect was added by hand. This happened because 
we did not use the ie prescription correctly: it had to be imposed in the propagators; this 
was already corrected in ref. ||10||. Taking into account the results from refs. fl and fllU|, we 



can conclude that the weak self-energy is null for an infinite stationary star. It should be 
repeated that following the correct ie prescription, when using propagators, is equivalent to 
taking the zero energy level for the states by matching to the free states far outside the star. 
Of course the definition of the zero energy level is essential when using eq. (§), as we shall 
now do. 

The null result for the (3 + l)-dimensional infinite star is a consequence of neglecting 
the existence of border effects. Now, before focusing on the spherical problem and in order 
to estimate the border effect, we concentrate on the geometrically easier problem: matter- 
free and neutron vacua separated by a flat border 0. For simplicity, the plane z = is 
assumed to be the flat border. The application of the approach resulting in eq. (H) requires a 
complete set of eigenf unctions for the Hamiltonian with the flat border. These eigenfunctions 
should be normalized following the same criteria as previously presented in section [TT]: they 
asymptotically behave as plane waves far outside the star, in the sense that they provide, far 
outside the star, the same probability and energy density as the free plane waves (solutions 
without a star): w(z = — oo) = 0. This condition is imposed inside each Hilbert subspace 



This problem can be understood as the consequence of zooming on the spherical border, the 
neutrinos wavelength being much smaller than the radius of the star. 
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with energy E G [E, E + dE}. 

Following [16|], we obtain the following set of negative energy eigenfunctions: 

^ e _ e (x) = N n _ \{e ik f x v h (k<) + e ik ?*Kv h (k<)) 9(-z) 
+ e lk ? x (l+H)v h (k<) 6(z) 
= N n _ [e^-(l + .7H(fe)0(-*) 



+ (e^ifcCfo ) + e^V« fc (^)) , (16) 



where they are directly related to incoming wave packets, the first coming from the free vac- 
uum and the second one from the matter vacuum; n_ labels the negative-energy eigenstate, 
n± = (±\E\, k p , h = —1); k p stands for the projection of the momentum on the flat border, 
and h is the negative helicity of standard neutrinos; iV n is the appropriate normalization 
factor. We also have fc > = (|i£|,Ar), k< = ( | _E | , A: < ) , and 

kf = (k x , k y , kf) , kf = (k X: k y , —kf) , (17) 
= (kx, k y , (kf)*) , fc^j = (k x , k y , —k^) ; (18) 

the z-components of the momentum are defined as 



k< = +\/E* - kl 



K = ^(\E\-W-kl- (19) 
The reflection coefficients JZ and J can be computed following, for instance, the work of 



Gavela et al. [16], we find: 1Z = — J = R^akp, where k p is the unitary vector in the direction 



of the transverse momentum k p , and a = 7°7; Ro is given by: 

2|fc|fe| 

fl ° " (*< + *>)'-»■ ' (20) 

If we consider a neutrino outside the star, with four-momentum (|£'|,/c < ), it can easily 
be seen that, for \E\ — \b\ < \k p \ < \E\, the matching condition generates a momentum 
inside the star with an imaginary component, k> = +i\Jk^ — (\E\ — |&|) 2 . Thus, a plane 
wave far outside the star in the above-mentioned parameter range generates a damped wave 
inside. We have, in other words, a non-penetrating wave. This is a crucial fact, because the 
probability density of waves inside and outside the star will be modified in a different way. 
It will be seen that the non-zero result for the weak self-energy, in both the flat and the 
spherical border cases, is dominated by the effect of these non-penetrating waves. 

Following Gavela et al., we can be sure that the only negative-energy non-penetrating 
solution, once n is fixed, is the first eigenstate of eqs. flltf ). Nevertheless, when transmission 
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occurs, both eigenstates are needed to span the whole eigenspace. It can be seen that 
eq. flT6[) appears to be a complete set of eigenfunctions [[U|, [jlTH . 

The normalization constants are chosen such that the probability density is asymptoti- 
cally equal to the free one. For a given energy this can only be achieved on average up to an 
oscillating term that will be considered as a negligible local fluctuation. This normalization 
convention will reach its unambiguous meaning in the next section. There, the asymptotic 
density of the solutions far outside the star is tuned to the one in the free vacuum. This 
asymptotic region is taken to fix the normalization because it extends to infinity while the 
star occupies a limited region. Applying then eq. we obtain the following weak energy 
density per energy: 

W\ {z) = -Wr ]E ^ riE ' z) - r ° {E ' z) } ■ (21) 

where T°(E,z) = 1, for free plane waves; T(E,z), for instance in the region \E\ > \b\, can 
be written as: 

T(E,z)=9(-z) + j(E)e(z) , (22) 

where 



-(E) = I <lr~== - ^ „ . (23) 



o VI -x 2 l &- 



The parameter b* is the dimensionless quantity b/\E\, and Rq is the same function Rq as 
given by eq. (|20D , expressed in terms of the variable x = \k p \/\E\ and of the parameter b*. 
It is easy to see that Ro is a real quantity, except for the non-penetrating waves, the result 
\Ro\ 2 = 1 being satisfied in the latter case. It is important to insist on the fact that the 
local fluctuations originated by the interference of incident and reflected waves have been 
neglected in both regions, z > and z < 0. The fact that these fluctuations contribute 
in a negligible way can be easily seen by performing the appropriate normalization of the 
eigenstates in a certain box. The contribution of the damped waves inside the star has also 
been neglected. That is why we integrate only for penetrating waves to obtain the function 
■y(E), i.e. the reason of the upper bound in eq. (0), which implies that Ro is always real in 
the l.h.s. of eq. (|23|) . The latter damped waves occupy a thin layer inside the star and are 
necessary to make the probability density on the border continuous. A similar effect will be 
discussed in the following section for the spherical case. 

In order to show the dominance of the effect of the non-penetrating waves, which is 



the main aim of this planar computation, the weak energy density ( f2~T| ) will be compared 



with a naive estimate of the one coming from non-penetrating waves: integrating simply 
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b|/E 
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FIG. 2. Plots of'y(E) (dashed line) and 7 (E) (solid line), defined in the text, as a function of 
\b\/\E\. 

the plane wave density in the two regions summed over all the allowed momenta in both 
k z = \J E 2 — \k p \ 2 for z < and k z = sj (\E\ — |6|) 2 — \k p \ 2 for z < 0, we obtain 



\E 




\h\ 




E 





(24) 



Equation (|24j ) results from computing separately the density of states in the star and in 
the free vacua, the eigenstates in these vacua being obtained for translationally invariant 
Hamiltonians, which do not generate reflection on the border. The difference between , y s (E) 
and 'j(E) is precisely due to this neglecting of reflection. The result fl24"|) only accounts for 
the fact that certain eigenstates, being allowed in the free vacuum, are not inside the star. 
There is an obvious one-to-one correspondence between these eigenstates and the ones for 
the untranslational-invariant Hamiltonian ([11]), which we called non-penetrating waves. 
In fig. |2|, the functions ^{E) and , y s {E) have been plotted. 

The small difference between the two curves in fig. |2] shows that the dominant effect of 
the border comes from the non-penetrating waves. This result will be confirmed in the next 
section. A few comments are in order here. First, the new feature in the (3 + l)-dimensional 
case with respect to the (1 + 1) one is that the non-trivial neutrino refraction index not 
only bends the penetrating waves and induces a reflection, but also, beyond the limiting angle, 
induces total reflection. This well-known phenomenon, acting here on the negative energy 
states, induces the dominant contribution to the star weak self-energy density. It must be 
stressed that this effect of non-penetrating waves is utterly unrelated to any Pauli blocking 
effect from the neutrino condensate inside the star. Indeed, as we have shown in the (1 + 1) 
case, in which the condensate exists, its effect is precisely to equate the energy density inside 
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the star to the one outside. On the other hand, the non-penetrating waves are reflected, not 
because of some states, which are already occupied, but because they tend to occupy states 
that simply do not exist, with an imaginary momentum. Finally, let us insist that 
eqs. ( PH ) and ( p3| ) imply that the border effect discovered here is a volume effect, i.e. it 
affects the energy density by an almost constant amount in the whole volume occupied by 
the star. This result is rather unexpected, a border effect being thought to act only on the 
surface. Of course it relies on the hypothesis that the wave packets extend coherently over 
the whole star. As will be discussed in section |VT[ such a hypothesis is perfectly sound for 
low-energy neutrinos. 



V. REALISTIC STAR: SPHERICAL 3-D 

We study a massless neutrino in the presence of an external symmetric static electroweak 
potential of finite range due to the interaction with the neutrons of the star. In order 
to calculate the weak self-energy J d 3 r w(r) = J d 3 r (<0|7Y(0)|6> - (0|ft (0)|0)), which is 
nothing else than the difference (H) integrated over space, we need to use the spherical 



Bessel functions [|18[ as a basis for the solutions of the Dirac equation. The effect of the 
static neutrons is summarized in the spherically symmetric square-well potential of depth 
b< (6 20 eV): 

.J b < for r < R 
^ ~ | for r > R ^ ' 

R being the radius of the star. 
^From the effective Lagrangian of eq. 

C eS = ^(i$-V{r) l0 ) V(r) , (26) 

the Dirac equation is: 

ce-pxp(r) = (E- V(r))ip(r) , (27) 

where 




Turning the kinetic energy operator a ■ p into spherical polar coordinates, we obtain the 
eigensolutions : 

^^("■fj™), (28) 
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where N is a normalization factor and X±k{@i 0) are the two-spinors |19] written in terms of 
the spherical harmonics Y™ and spin 1/2 eigensolutions xi 



X K >o( 



53 («a*|m. («-i)m)y^,. 



/ 



2/t+l 1 K,m— § 



2k+1 K,m+2- 



(29) 



x™<o( 



-ft 





l 


— ft+m 


2 


-2k- 


1 


— k— m 


1 
2 


-2k- 


1 







(30) 



re— l,m+i y 



Here X±| K |(0>0) are eigensolutions of J 2 , J 2 , L 2 and S 2 , where J = L + S, J z , L and 
S are respectively the total angular momentum, its projection along the z-axis, the orbital 
angular momentum and the spin angular momentum verifying the eigenvalue equations: 

(*.L + 1)X£ W (M) = =FNx£|„|(M) , 

J, X £ W (M) = ™X? N (M) , -J<m<+J, (31) 
|/c| = (1,2,...) , |/e| = J+l/2. 

With this notation, we can easily check that the two chirality states are described by the 
same spinor of eq. (|28|): 

1 



1>Z.n = 2 (1 T 75 ^ 
the positive and the negative values of k being related through the relation: 

• _, „i,m „l,m 
I 75 W-\k\ = V\k\ ■ 



(32) 



(33) 



It then suffices to compute ipT f° r positive values of k, use eq. ([33|) for the negative 
values of k and then get the left- chirality neutrino state using eq. (|32"D . Let us from now on 
write the positive k's as k = I = 1, 2, 3, with J = I — 1/2 (see eq. (|29|)). 

The radial functions Uxp defined in eq. (p8|) verify the following coupled differential 
equations: 



d 



r uur) 



dr 1 r «2(^) 



V (E-V(r)) 
<E - V(r)) I 




(34) 



Now, decoupling the equations and using eq. (|25f) and the fact that the solution must 
be regular in the centre of the core, i.e. r — > 0, we get: 



Ul (r) = a [ ji (k r) 6(R - r) + (B 3l (E r) + C Vl {E r)) 9{r - R)\ 
Mr) = a[ (k r) 6{R - r) + (B j^E r) + C y z _ x (£ r)) 6{r - R)] 



(35) 
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where we define 



k = E - b, 

a = = ; (36) 

VB 2 + C 2 V 7 

B, C being integration coefficients to be fixed from the matching condition at the surface. 
The factor a will be justified below to guarantee the good asymptotic behaviour of functions 
Ui^ij)- The functions ji and yi are the spherical Bessel functions of the first and second 
kind. In the inner region, the solution can be written only in terms of ji because the wave 
functions have to be regular at the origin and the j/z(p)'s are not (see eq. flSOP). 

According to eq. (||), we are going to use only negative energy states and thus negative 
arguments for these Bessel functions (E, k < 0) and the following relations are helpful: 

m = (-i) m jK-p) , vi{p) = (-i)'w(-p) • (37) 



The solutions Uip have been written in the form fl35|) in order to match, up to a phase, 
the free solutions as r — > oo, where the asymptotic behaviour of the Bessel functions is: 




ji(p) p^oo sin \p- — J jp , y t (p) p^oo cos - — ] jp . (38) 
This behaviour implies the one followed by our solutions U12' 

Ul {E r) r ~ D 




/ \ / 7rZ 

u 2 [E r) r J^oo cos I E r - — 
where <fi = arctan(C/5) . (39) 

Needless to say that in the absence of the star (b = 0), a = 1 and the spherical solutions fit 
exactly the free ones, in the asymptotic region r — ► oo. 

Both inner and outer solutions must join at r = R, and this matching fixes the coefficients 
B and C as follows: 



B 



C 



j t (k R) y^jE R)-Ji^(k R) yj (E R) 
j t (E R) y l _ 1 {E R^-j^iE R) m (E R) ' 

(40) 

j;_i(fc R) ME R)-ji-i(E R) gk R) 
j t (E R) yi ^(E R)-3i-x{E R) y t (E R) ' 



Finally, in order to fix the solutions definitively, we need to calculate the normalization 
factor iV by following the asymptotic normalization convention already introduced in the 
previous section: we impose the probability density to be asymptotically equal to the free 



17 



one. Since we have fitted a to make the solutions with b 7^ coincide, up to a phase, with 
those of the free case (b = 0) when r — > oo, we need to compute N only for the free case. 
In the latter case, b = 0, a = 1, the Hilbert space of spherical solutions is the one of the free 
plane waves j^s^' r - The solutions (see eq. (PH)) reduce to 



N 



j t (E r) xT{9A) 
iji-^E r) x™t(9,, 



(41) 



Now we compute the density of these states p(E) in the Hilbert subspace corresponding to 
E G [E,E + dE], by summing up all the quantum numbers and the angular degrees of 
freedom, while keeping the left chirality and spin fixed because we are counting the negative- 



energy and left-handed neutrino states. This gives the probability p(E)dE 



(2tt) 3 



the relation 18 



^(2/i + l)|j M | 2 = l, 

jtt=0 



dE. Using 



(42) 



we obtain: 



N = E/V2^. (43) 
Now that the solutions of the Dirac equation are known, we can calculate numerically 



w(r) } which can be written following eq. @ as [14]: 

w(r) = f dE J2 E (Vf t(*0 Vf(*0 - € (0) t(r) ik { °\r, 



(44) 



Lm 



where the sum is over all the degrees of freedom of the negative-energy left-handed neutrinos 
and where the subscript (0) refers to the matter-free solutions (6 = 0). The integration over 
the angular part and the summation over m is direct, using eqs. (|29| ) and flSOP , which lead 
to 

l r ,s ,s (2J+ 1) _ 21 

Air ' 



Ait 



dnj2x\r(e,m' 



Att 



(45) 



The weak self-energy density (w(r) = dW/dV = 4 ^ dr , W being the total weak energy and 
V the volume) is then simplified to (using the spherical symmetry); 



w(r) = f° dE^Yl KPi(k, r) - pf ] (E, r)) 



(46) 



where we have defined 

pi(k, r) = a 2 { [jf(k r) + jf^k r)] d(R-r) + 

[(Bj^E r) + C m (E r)f + (Bj^E r) + Cy^E r)) 2 ] 6(r - R)} , (47) 
p\ 0) (E,r)=jf(E r)+jf_ 1 (E r) , 
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k and a being given by eq. fl36|). We were not able to perform the summation over all 
the values I and negative energies analytically. We then use numerical computations of the 
Bessel functions, and use some knowledge (see ref. |TB[) about these functions in order to 
justify the truncation of these infinite summations in I, as will be explained now. 
Let us define the function 

f(E,rJ 2 ) = j2KPi(k,r)-pf >) (E,r)) , (48) 
1=1 

We have observed that for r > R and fixed values of E and the result becomes rather 
stable from r > R up to r 3> R, while there is a sudden change when we pass the border 




r = R. When we integrate the density over r, the main contribution to (|48|) comes from 
inside the star, re [0, R] and thus the effect is grossly proportional to the volume. 

In order to understand these results qualitatively, a few remarks about the behaviour 
of Bessel functions are appropriate. The spherical Bessel functions are solutions of the 
differential equation: 

-M)(*>-, (J + *-ffi±I>)<„>-o. W 

For p large, say p ~ po, we may get a hint by replacing 1/p 2 by 1/p 2 , m eq. fl49f) . 
For A; 2 — /(/ + l)/po < 0, i.e. Z <; |&| po, the solution j; (y;) is a damped (exploding) exponen- 
tial []. This damping (explosion) is expressed by the well-known behaviour of the spherical 
Bessel functions ji and yi in the vicinity of the centre: 

*M-°(2fTiw ■ ^^fsrrr ( ' +,) • (50) 

For k 2 — Z(/ + 1)/pq > 0, - i.e. I ^ \k\ po, the solution is an oscillating function as ex- 
pressed by the asymptotic behaviour (|38|). 

In practice the transition between these two asymptotic regimes is rather fast. In other 
words, the Bessel function j t (p) increases with p at first as a power of p, then increases 



exponentially until a transition at p = 1) where it becomes an oscillating function, 

see eq. (|38|). Spherical Bessel functions are not so different from 



4 This damping of ji expresses the semi-classical fact that I ~ \\po x k\\ < po\k\, the waves corre- 
sponding to I > po \k\ are damped. 
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j7(p) = Hp - y/Ki + 1)) sin ( P - j ] / p, v P , 

Ztt 




yr(p) = % - V / ( / + !)) cos ( p - y ) /p. p- V*0 + 1) > o. (r,i ) 

The exploding part of yi is not described by eq. (0), but we do not care since, precisely 
from imposing regularity at the origin, the yi's only contribute to our solutions outside the 
star, eq. fl35|), where they are in their oscillating regime, as can be easily checked. 

In the following we will use this approximation to guess the results and we will check 
numerically these guesses and estimate the corrections to these approximations. Such a 
strategy is necessary since we are totally unable to calculate numerically with the realistic 
value of \b\R, which is of the order of ~ 10 12 . This would need computing up to angular 
momenta larger than 10 12 ! In practice we have modestly limited ourselves to \b\R — 10 
and \E\/b going up to ~ 100. It is then mandatory to get a qualitative understanding of 
the results to be able to extrapolate them to the realistic values. Before going further, it is 
worth remarking that the function defined in ( |5l|) verifies the saturation relation fl42|) at 
leading order in l/(|fc|r). 

From ( fTi]) it seems that the summations on I may be truncated at I > \E\r. In order to 
check this conjecture we will use the free case in which we know the exact solution. For the 
plane waves, the energy density w p due to left-handed negative-energy neutrinos is 

w p {r) =JdE — = JdE—Y, l P f ] (E, r) , (52) 

as can be directly computed from the plane wave functions, and also via the definition of 
p[ ^ in (f47|) with the help of relation (0). The integration over E in eqs. (|46| ) and (|52|) is of 
course divergent. We will come to this ultraviolet problem later, in section [VI]. 
We have computed 

dG p f R „ o , E 3 

= / 4vrr 2 dr — - 53 
dE Jo 2tt 2 v ; 

and compared it with 



?3 h 



dG C R . o , E 

i=i 

in order to see the convergence of the summation. It is confirmed that the summation 
Y^i=xlp\ i(E, r) is already very close to its limit 1 as soon as l 2 > \E\ R. This can be seen in 
fig. 3, where we have plotted the absolute values of eqs. ( [5*3"D and (|5l|) . 

We have also checked that the summation in (|48|) saturates when \E\ r < l 2 . 

The numerical computations show that a, eq. (0), is quasi-independent of I for 
I < \E — b\ R and is verified for: 
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FIG. 3. 



in eq. n5l) as a function of the energy expressed in units of \b\. The summation 
on I goes from 1 to the value indicated on the figure, r is integrated between and R (the radius R 



dG 
dE 



of the star has been fixed to 10|6 The bold line represents the plane wave density 



dG z 



dE 



in eq. 



a 2 \E-b\R>i a 2 t = y E j • (55) 

This result may be understood by the following argument. Using the approximation 
(|5"T|), the waves that do not vanish inside the star near the border have the form of a sine 
divided by (k R). The corresponding ones, outside the star, are a combination of a sine and 
a cosine divided by (E R). Taking the trigonometric functions to be of the same size on 
average, an obvious factor of a ~ a t = \k/E\ is needed for the matching at R of the inside 
waves and the outside ones. Interestingly enough, this result leads, combined with (^2]), i.e. 
summing over all angular momenta, to an average energy density inside the star of 

E(E - 6)7(2tt 2 ) ' (56) 

which is the energy multiplying the plane wave probability density, see eq. (|52|), but for a 
plane wave shifted by the potential b. This is reminiscent of eq. (p4j). 

The validity of eq. (|5^) is shown in fig. [|, where we compare a to a t for different values 
of the energy and the parameter /. We can see on this figure that a 2 — a 2 is negligible in 
average, as long as the momentum / verifies I <C \E — b\ R. 

At this point, we now have a "naive" or "crude" formula for w(r): 
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FIG. 4. a 2 — a 2 (see eqs. (51 ) and (jgqjj as a function of I, for different values of the energy. 
The star radius is R = 106 -1 . 



w crude (r) = JdE^((E~ b) 2 - E 2 ) 9{R - r) , (57) 

of course, this is a first approximation, the same as we had reached in the flat case (see 
section [TV]) from the function •y s (E), eq. (21). 



The reason for this agreement between the flat case and the spherical one is easy to 
understand physically. It is related to the non-penetrating waves. Let us first consider the 
waves with I < \E — b\ R. We have already argued that they are trigonometric functions 
divided by \E — b\ R (\E\ R) inside (outside) the star. They are matched on the border. 
From that matching itself and their oscillating behaviour, their density has to match up 
to a large distance from the border. A different and crucial effect comes from the angular 
momenta 

\(E-b)\ R< I < \E\ R . (58) 

In the trigonometric regime these waves are functions divided by \E\ R outside, while 
inside they are in the fast decreasing regime. The matching adjusts them for r = R, but, as 
r decreases inside the star, they fall very fast. Essentially they are non-penetrating waves 
exactly analogous to the ones with \E — b\ < k p < \E\ in the previous section. Indeed, 
since L ~ r x k, and noticing that the radial vector being perpendicular to the surface, 
its cross product with k is parallel to the surface, it becomes obvious that the inequality 
\E — b\ < k p < \E\ is equivalent to eq. Q58|) . 
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We have thus seen qualitatively, and checked numerically, that the dominant contribution 
to w(r) is due to non-penetrating waves in the spherical case as well as in the planar case. 



The numerical check consisted in the following steps. 
In the beginning, we compared: 



dW 
~d~E 



= / ^dr — Y,l{ P i{k,r)-pf\E,r)) 

JO Z7l l=1 



(59) 



for several values of I2, r max being greater than the star radius R, with 

dW cru de E 



dE 



2tp 



[(E-bf-E 2 ) 4vri? 3 /3 



(60) 



In fig. 5, we have plotted eq. (^Sj) as a function of the energy for several values of I2 
and compared with eq. (|60"D. The integration over r in eq. (|59| ) has been done from to 
r max — 3-R. This figure confirms that, as long as \E\ R < I2, there is no sizeable difference 
between the exact result fl59D and the "crude" one (|60|). For higher values of \E\ R, we must 
increase I2 in order to take all the contributions into account. 




-9 -8 -7 -6 -5 -4 -3 



-1 



FIG. 5. 



dW 



dE 



in eq. (\5Q) as a function of the energy expressed in units of\b\. The summation on 



I goes from 1 to the value indicated in the figure, r is integrated between and r n 



30 1 6" 1 1 (the 



radius of the star has been fixed to 10\b The bold line represents the energy density 



dW r . 



dE 



in 



eq. tfiQ). 



After that, in order to check the corrections to our formula eq. 
computed: 



more accurately, we 
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d5{r) _ d6in d5 out 
dE ~ dE dE 



where 

ddiJr) 



dE ~ 4?Tr 2 - 2 



E(a 2 -a 2 t )l(jf(k r)+jf_ 1 (k r)) 
1=1 



0(R-r) 



and 



(61) 



= Anr 2 ^ £ I [a 2 ((B 3l (Er) + C Vl {Er)f 



dE 



i=i 



+ (Bj^Er) + Cy^Er)) 2 ) - p^(E,r)\ 6(r - R) . 

Here, we are subtracting from the exact formula the "crude" one, expanded in Bessel func- 
tions with the help of the constant factor a t . Indeed, from the plane wave expansion ( |52|) 
and the definition of a 2 fl55l), it is obvious that the term proportional to a 2 in (|6TD sums up 
to the one proportional to EiE — b) 2 in flBTjf ) and that the term proportional to E 3 cancels 
the sum of the terms p\ ; . We have expanded the "crude" contribution into Bessel functions 
in order to minimize the error due to the truncation of the sum. 
After having performed the integration over r: 



dA 

dE ' 
where 

dAir 



dA iTi 
dE 



+ 



dA 



out 



dE jo 
we compare in fig. || eqs. ( 



dr 



dE 

dd in 
dE 

and (|6 



and 



dA 



out 



dr 



d5, 



(62) 



out 



dE Jr dE 
We can see from the figure that the correction to 



our naive formula ( pTQ is relatively very small and better for larger energies 
The integration over r, where r max > R, shows that 



dA, 



< and 



dA 



out 



> o, 



dE dE ' (63) 

which means that our crude estimate is an underestimate (overestimate) of the probability 
density inside (outside) the star. This fact is due to the tendency of the exact solution 
to provide a continuous matching on the border, while our "crude" formula has a jump. 
We shall discuss this in more detail. Moreover, we have found that this overestimate and 
underestimate approximately cancel as can be seen from fig. |^. Note that in this figure 
(fig. 0), the summation over I goes to I2 = 80, the fluctuations that we see for E < — 8|6| 
are due to the truncation in the sum over I when \E\ R > 80. 

Let us look in more detail into how the sudden jump of our naive formula ( |57|) when 
crossing the border of the star is smoothed down into the exact result. There must be some 



'The falling down of the dashed line in fig. |6| for E < — 8|6| is a truncation artefact. 
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FIG. 6. The dashed line represents 



dW 



dE 



in eq. fibl\) and the solid line 44 in eq. (\63l). They 
are functions of the energy expressed in units of \b\. The summation over I was done from 1 to 80 
and the integration over r from r = to r = 30 1 1> 1 1 (the star radius being R = 10 1 £> 1 1 ^ - 



layer in which the transition takes place. As we have already claimed, this layer is dominated 



by the functions in the domain of eq. (|58|). The width of the border effect depends directly 
on the speed of the falling-down of the Bessel functions in (^). If we use the fact that, for I 
high enough, the spherical Bessel functions behave like (\k\ r) 1 (see eq. (pU"D), we can try an 
estimate of the border effect, which gives the correction to our "naive" eq. (|57|) . The mean 
value of I inside the domain (f)8|) is 



R 




(64) 



If we suppose that half of the jump occurs at r = _R_ < R in the interior of the star and 
the other half outside r = R + > R, we can construct a term proportional to 



dS{r) 
dE 



E*_ 
2tt 2 



rr> \ 2/ mean 

r) 



9{R-r)- l -{l 



a. 



-)' 



6(r - R) 



(65) 



which will have the role of smoothing the jump in our naive approximation formula eq. 

In order to confirm this effect, we have compared in fig. |8], for different values of the 
energy, dS/dE with d5/dE. 

In fig. ||, we can see the agreement between the l.h.s. of eqs. (|65|) and (6lf), except for 
E such that \E\ R> l 2 , where the truncation effects are important. 
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FIG. 7. The horizontal axis represents the energy in units of \b\. The summation over I has 
been performed from 1 to 80. For E > 8b, the positive curve represents (eq. the contributions 

dA ou t I dE of the exterior of the star, the negative curve the contributions of the interior. For 
the exterior contribution, r have been integrated from R to 3R, the star radius R being fixed to 
R = 10 1 1> 1 1 . An accurate cancellation of both pieces is seen. The oscillation seen for E < 8b is 
due to the truncation artefact. 



To have an estimate of the width d in which the "joining" occurs, it suffices to write 



ry \ ^mcan / <y 



u j exp \ l meaa !og(^ 

defining d such that the probability decreases by a factor of 1/e, we get 

d= whm- (66) 

From fig. we can easily see that the larger \E\ is, the narrower the width of the joining 
layer is. This is clearly seen in fig. |8| a while in |] b, beyond E = 8b, the truncation effect 
comes into play, since \E\ R > l 2 — 80. 
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dS 



eq. jjfiBQ and the solid line in eq. (61), as 



functions of the position r, in units of (the radius of the star has been fixed to R = 10 \b~ \) 
for different values of the energy. The summation over I has been done from 1 to 80. The two 
lines are similar, confirming the validity of approximation ffih\), except for E = 10|6| where the 
truncation artefact dominates. 



VI. DISCUSSION 

We now have a qualitative understanding of the dominant contribution to the star mass 
correction due to the neutrino exchange. It doesn't vanish when we consider the realistic 
(3 + l)-dimensional star, because of a border effect proportional to the volume of the 
star. It is ultraviolet- divergent and we have to get some understanding of the ultraviolet 
cut-off C. This cut-off corresponds to the limit up to which our theory describes, to a good 
approximation, the exact one. The dominant contribution (see eq. (fT7D) w cru d e (r) to the 
total weak energy is equal to (see eqs. (|_ 

4 



W crude = ^R 3 f° dE E((E - b) 
Sir J-c 



and (0)): 
E 2 ) 



9?r~ 



-R 6 C% 



(67) 



This result presents the unexpected feature to be odd in b, while in the perturbative cal- 
culation the contributions with an odd number of neutrons vanish exactly, resulting in a 
total energy even p] in b . It is worth noticing that had we started from the symmetrized 



3 We are indebted to Ken Kiers and Michel Tytgat for having pointed out this fact to us. 
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expression for the vacuum loop eq. @, we would have obtained, instead of (|67j). 

W sym _ crude = ^-R 3 J° c dE 2Eb 2 ~ ^R 3 C 2 b 2 . (68) 

The effective Lagrangian ([|) is valid only under the assumption that the neutrons are 
approximately static. This breaks down, of course, for an energy scale of, say, C ~ 100 
MeV. Beyond that scale, the neutrons feel the recoil of the scattering. Not far above, one 
encounters the scale of confinement in QCD ~ 1 GeV, and the interaction start to "see" the 
substructure of quarks and gluons. 

At about the same scale, C ~ ^, where r c ~ 0.5 fm, as noticed by Fischbach a 
repulsive core prevents the neutron from "piling-up" in space. In the same ballpark, the 
average distance between neutrons 1/C = n~ x ' 3 ~ 1 fm [n n is the neutron density in 
the neutron star, n n ~ 0.4 fm -3 ) is the inverse of an ultraviolet cut-off since, for smaller 
distances, our picture of an homogeneous background of neutrons breaks down: the neutrino 
"see" individual neutrons in a vacuum that is the standard vacuum. All these cut-offs are 
of the same order of magnitude. Choosing the latter, C = n l J 3 ~ 1 fm -1 , we get from fl6"7|) 
per unit volume an energy of 

bC 3 \b\n 

W crude ~ ~ ~ 10~ 9 GeV/fm 3 ~ 10~ 8 m n per neutron, (69) 

3tc 20 

m n being the neutron mass, which means an interaction energy eight orders of magnitude 
below the neutron mass, and thus totally negligible. From fl68|) , the result would be even 
smaller: 

^2/^2 \b\ 2 n 

W sym - crude ~-—~- ] -^~-\Q- 1G GeV/{m 3 ~ -10- 15 m„ per neutron. (70) 

We did not yet manage to understand theoretically the relation between the results ( |69"|) 
and flTUp. They are derived respectively from eqs. (|8|) and (0). On the one hand, eq. ([|) 
leads to a result even in b, as does the perturbative expansion; on the other hand, taking the 
negative energies for the vacuum seems more physical to us. It should also be noted that eq. 
(^) is equivalent in Fourier space (see eq. (|T0|)), to averaging the complex integration closed 
in both the upper and in the lower half-planes. The issue is clearly related to the ultraviolet 
regularization. Indeed, subtracting (||) from (||) and using the closure theorem, it is easy to 
obtain the formal result: 

i d 

w(x) - w sym (x) = - -^rh(x - y) = o , (71) 

where the vanishing occurs since we derive on time a time-independent quantity. Once 
regularized by an ultraviolet cut-off so that — C < E < C as in ( |67D and (|^f), w(x) —w sym (x) 
is not zero. We therefore conclude that the difference between (|69| ) and ( |70|) is related to 
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different ultraviolet regularizations. At this point we prefer to present both results, while 
studying the issue further. 

Another argument leads to reduce the estimate (|oTS|). Indeed, throughout this paper 
we have considered stationary eigenstates of the Hamiltonian that extend over all space, 
such as plane waves. The use of such extended waves contains the implicit assumption that 
the neutrino wave packets extend over the whole star and far beyond. In other words our 
solutions know about the whole star and its surroundings. This is perfectly legitimate for 
low-energy neutrinos, since their mean free path is much larger than the radius of the star. 
This means that they feel the coherent interactions from the neutrons as expressed by the 
effective Lagrangian (^), but they do not experience incoherent scattering with the neutrons 
of the star. 

A rough estimate of the neutrino cross section with neutrons is a ~ E 2 G 2 Fl where Gp is 
the Fermi constant and E v the neutrino energy. The mean free path is approximately 

1 1 /0.5MeV\ 2 , Zl00keV\ 2 „ rl 
A ~ — — r ~ km ~ 25 km , 72 

n n a n n ElG\ \ E v ) \ E v ) y ' 

where n n is the neutron density in the star. This means that for an energy of 100 keV, 
the mean free path is ~ 25 km, of the order of the star radius. Let us call this auxiliary 
ultraviolet cut-off c (c ~ 100 keV). 

We now decompose the integral on the negative energy modes in two parts. For \E V \ < c 
we use formula flSTD - For energies larger than c but smaller than ~ 0.2 GeV, we may argue 
that the wave packets must have a size of the order of the mean free path. Therefore the 
"knowledge" of the border effects extends inside the star on a distance of the order of A 
from the border. A neutrino deeper in the star does not feel the border, we are back to the 
situation of an infinite star |4[] where we had a vanishing result. The non-zero contribution 
comes only from a thin region around the border of width A. Instead of R 3 in (|67D , we 
take R 2 X ~ (c/E v ) 2 R 3 (using the 1/E% law in ( |72|) and the fact that for c = \E U \, we have 
A ~ R). We thus have, instead of 



-R 3 I dE {-2)E 2 b +-R 3 I " dE{-2c 2 )b = -R 3 b{c 3 /3 + c 2 (C - c)) ~ R 3 c 2 Cb . (73) 

T J-c 7T J-C IX 

Since c/C ~ 10~ 3 , we get an average estimate for the energy per volume still six orders 



of magnitude below the result of (|69D , i.e. 

~ 10~ 14 m n per neutron. (74) 

Finally it should be remarked that all along this work we have solved a stationary prob- 
lem. This means that we have assumed the star, and also the neutrino states inside it, to 
have reached an equilibrium status. It is well-known that the stars evolve during their life; 
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we thus have implicitly assumed an adiabatic adjustment of the neutrino states. Since the 
star evolution is slow and since neutrino motion inside them has the velocity of light, this 
seems to us a reasonable assumption. Some further study might still be welcome. We have 
also assumed a zero temperature for the neutrinos, since we believe that their interaction is 
too small for them to thermalize. 



VII. CONCLUDING REMARKS 

In order to settle definitively the question of the stability of a neutron star, the multibody 
exchange of massless neutrinos has been computed analytically and numerically for a finite 
star. 

The effect of a border is twofold. First it induces in a natural way the neutrino condensate 
as proved in The latter condensate, does not produce any neutrino exchange interaction 



energy in the simplified (1 + l)-dimensional case and we find it to be negligible in the realistic 
(3 + l)-dimensional case. 

The second effect of the border is that the neutrino zero point energy inside the star differs 
from the outer one because of negative-energy waves that cannot penetrate inside the star, 
being beyond the limiting refraction index at the border. This contribution is proportional 
to the volume of the star, but it is still tiny (10~ 8 -10~ 13 GeV per neutron), completely 
negligible in comparison with the neutron mass. We find no infrared divergences in the full 
non-perturbative result, which would have necessitated the introduction of a neutrino mass. 

The general conclusion of this work is that the neutrino does not need to be massive to 
ensure the stability of a neutron star. This is in agreement with recent works (commented 
below) by Kachelriess and by Kiers & Tytgat ||. There is no catastrophic effect due 
to the multibody massless neutrino exchange. As already stated in refs. [|J and ||10|| , this 



catastrophic result claimed in [|7| is only due to an attempt to sum up the perturbative series 
outside its radius of convergence. 

While finishing this paper, there appeared a paper by Kiers and Tytgat [fj. They accept 
the point of view developed in Q for an infinite star and try, as we do, to solve the problem 
of a finite star. Their starting goal is, as ours, to compute the density in eq. (|8|). They 
use a clever technique based on quantizing in a large sphere and expressing the vacuum 
energy density in terms of the phase shifts. They first study analytically the unphysical but 
illustrative case of small bR and then numerically the large bR case. They show that the 
perturbative series a la Fischbach already breaks down as early as bR > 7r (for the neutron 
star, bR ~ 10 12 ) while the non-perturbative calculation gives a negligible result, a conclusion 
which we fully share. One difference between their result and ours is that they find a relation 
between the energy density of the neutron star and that of the neutrino condensate. We 
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find on the contrary that the result is mainly due to non-penetrating negative-energy waves, 
which are not related to the condensate. We did not yet manage to understand the reason 
for this discrepancy. It might be related to different UV-regularization methods. 

Kachelriess [H also agrees with us about Fischbach's "catastrophic" result. He computed 
the total weak self-energy for an infinite neutron star following the Schwinger method, by 
using the neutrino propagator in momentum space, as we did in ref. H. He obtained 
a non-zero weak self-energy without taking into account the neutrino sea effects. As we 
acknowledged in section |V| and in ref. flO ], a minor mistake was made in ref. Q: the 
contribution of a pole had been forgotten in the calculation. Once this mistake is corrected, 
we agree with Kachelriess about the result when neglecting the neutrino sea. He attributed 
the discrepancy to the fact that we took the limit (y — > x) before integrating over the whole 
space to obtain the total weak self-energy (see eq. (Q)). We do not agree with this conjecture 
about the discrepancy, first because, as we just mentioned, the forgotten pole removes the 
discrepancy, second because we have verified that our previous UV-regularization makes 
the x — > y limit regular. Finally, from our analysis of finite stars we insist that the 
condensate has to be taken into account and, amazingly, it exactly cancels the forgotten pole 
contribution, resulting finally in w(x) = for an infinite star. 

A few weeks later, Fischbach and Woodahl [|20| repeated Fischbach's original claim and 
used the same expansion, order by order, in the number of neutrons, i.e. equivalently in 
perturbation in the parameter bR. Astonishingly enough they did not consider the series 
of works demonstrating that this series is simply divergent, but that the total result may 
be computed directly by the effective Lagrangian technique. They argued that our non- 
perturbative calculation encounters cancellations because, in the effective approach, the 
neutron medium is assumed by us to be a continuous background. Of course, treating the 
neutron medium as a homogeneous continuum medium is an approximation a la Hartree- 
Fock, and it should be corrected in the ultraviolet by taking into account the correlation 
between neutrons. This is precisely one of the reasons why we considered that a natural 
ultraviolet cut-off was the energy scale of a few MeV. The authors of take the size of the 
neutron hard core as an ultraviolet cut-off. Why not, although many other ultraviolet effects 
arise at the same scale of a few 100 MeV: the neutron recoil, the quark and gluon content 
of the neutron, without forgetting the incoherent neutrino-neutron scattering discussed in 
the previous section. But these ultraviolet effects will not at all modify the analysis of the 
infrared catastrophe advocated by Fischbach and denied by us. The authors of 0] seem 
to imply that we have added some unjustified assumption in our work. The truth is on 
the contrary that we have assumed nothing that they did not assume themselves, such as 
the static neutron assumption, but we have not assumed, as they do, that a neutron is not 
allowed to interact more than once, neither did we make the drastic approximations that 
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appear in their work at high order in perturbation. The effective Lagrangian approach allows 
to compute exactly, in a simple manner, and with fewer assumptions than the perturbative 
expansion approach. 

The authors of [E0| criticize our recent (1 + l)-dimensional toy calculation [TO] arguing 



that the critical parameter bR is much smaller than 1 in (1 + 1) dimensions. However, they 
did not notice that our (1 + l)-dimensional result is absolutely exact, independently of the 
parameter bR, which, incidentally, we have taken to be large. 

To finish, we feel it necessary to insist. The main issue is the failure of the perturbative 
expansion, which is infrared- divergent. Happily one can spare this difficulty thanks to the 
effective action technique. Once this point is understood, the different analyses all agree, 
notwithstanding minor discrepancies, that although the massless neutrino exchange between 
fermions is a long-range interaction, it does not give any significant contribution to the total 
energy of a neutron star, finite or infinite. 
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